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Abstract. We report the finding of two large-scale structures of galaxies in a 40 x 35 arcmin 2 field embedded in a 25 deg 2 area 
where two 100/r'Mpc-scale structures of quasars and Mgn absorbers at z — 0.8 and z — 1.2 have been previously detected. 
Using deep optical (V- and /-band) imaging, we are able to select high-redshift (0.5 <z< 1.3) early-type galaxies as those 
redder than the cluster red sequence at z = 0.5 and having / < 23.5. Through comparison with a 35 x 35 arcmin 2 control field, 
we find a 30% excess of these red galaxies, corresponding to 563 galaxies across the field. The colour distribution of the excess 
galaxies shows a coherent peak at 2.7 < V - / < 3.1 and the magnitude distribution is well described by a Schechter function, 
showing unambiguously that the galaxy excess is primarily due to a population of early-type galaxies at z = 0.83 ± 0.08. In 
follow-up NIR imaging of the four best cluster candidates, we identify three clusters at z = 0.8 ± 0.1, each having well defined 
red sequences of ~ 20 galaxies at / - K s — 3.4, forming a structure 20/r'Mpc across. The fourth candidate, corresponding to 
a previously identified cluster at z = 1.2 ± 0.1, is confirmed with 28 K s < 20.5 galaxies having VIJK S colours consistent with 
passively-evolving galaxies at z - 1.2 located within 1 arcmin of the cluster centre. In all four NIR fields we find numerous 
galaxies with the VIJK S colours of passively-evolving galaxies both at z — 0.8 and z - 1.2, including regions well away from 
the clusters. The results of this study support the presence of two sheet-like large-scale galaxy structures extending across the 
20-30 /r'Mpc scales explored by the optical data, one "sheet" at z — 0.8 with three associated clusters, and a second "sheet" at 
Z—1.2 with one embedded cluster. Our results confirm that, at least inside the probed region, the two known 100/?~'Mpc-scale 
structures at z ~ 1 of quasars and Mgn absorbers also mark out volumes with an enhanced density of galaxies. The finding of 
two galaxy structures at z — 1 on scales of 20-30 /?~'Mpc within much larger superstructures of quasars and Mgn absorbers, 
firstly establishes the presence of galaxy sheet-like structures at z ~ 1, and secondly, suggests that they may extend up to 
100/^'Mpc-scales. 

Key words. Cosmology: large-scale structure of the Universe; - Galaxies: clusters: general; Galaxies: evolution; - Quasars: 
general 



1. INTRODUCTION 

Over the last twenty years galaxy redshift surveys have suc- 
cessfully mapped the local universe, revealing that on scales 
of > lOO/z^Mpc the spatial distribution of galaxies appears 
to be cellular or sponge-like, with patterns of wall-like super- 
structures having typical extents of 30-50 /; _1 Mpc and thick- 
nesses of ~ 5 h~ l Mpc, surrounding low-density regions or 
voids 50-70 /i _1 Mpc across (Doroshkevich et al. 1999, 2000, 
2001). Recent surveys covering representative volumes of the 
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universe have established that the sheet-like large-scale struc- 
tures such as the "Great Wall" of galaxies (Geller & Huchra 
1989), that were the most striking features of early surveys are 
typical phenomena (e.g. Shectman et al. 1996; Doroshkevich et 
al. 1996, 2000, 2002), with average extents of 50-100 /z^Mpc, 
comprising ^ 50% of the galaxies with an overdensity of ~5- 
10 above the mean (Doroshkevich et al. 2000), their distri- 
bution within these sheets being inhomogeneous with galax- 
ies concentrated in clusters and filaments (e.g. Ramella et al. 
1992). There are, nevertheless, indications for deviations out 
to scales of ~ 160/i _l Mpc (Best 2000), and a turnover in the 
galaxy (cluster) power spectrum at these scales is still disputed 
(Miller & Batuski 2001). The "Great Wall" appears to be rep- 
resentative of the upper limit in size (~ 100/i _1 Mpc) for these 
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sheet-like entities belonging to the observed super large-scale 
structure, and as such they are expected to be rare. 

The "standard" model for the origin and evolution of the 
universe is the Cold Dark Matter (CDM) model (Blumenthal et 
al. 1984), whose most recent version has passed stringent con- 
sistency tests provided by analysis of the CMB, supernovae, 
and galaxy surveys (Spergel et al. 2003). The ACDM model 
together with observational constraints implies that ~ 90% of 
the mass is dark, thus the distribution of luminous matter such 
as that of galaxies is superposed upon a dominant background 
DM distribution. Our ignorance of the relation between spe- 
cific luminous markers (i.e. intergalactic gas, galaxies, radio 
galaxies, quasars) and the underlying DM distribution is ex- 
pressed through a "bias" parameter or function (Kaiser 1984; 
Blanton et al. 1999), which represents one of the most funda- 
mental sources of uncertainty in the efforts to perform mean- 
ingful comparisons between the observations of the galaxy dis- 
tribution and the DM model predictions. 

With these caveats in mind, the observed large-scale struc- 
ture in the galaxy distribution may be understood in the context 
of the Zel'dovich non-linear theory of gravitational instability, 
whereby perturbations in the initial Gaussian random density 
field collapse rapidly along one axis to form Zel'dovich pan- 
cakes (sheets), with further collapses along the remaining axes 
to form filaments and clusters (Demianski & Doroshkevich 
1999a, 1999b). Numerical simulations have allowed the evo- 
lution of large-scale structure from the initial Gaussian per- 
turbations to be followed up to the present for a variety of 
CDM models, successfully reproducing the main properties of 
the observed large-scale galaxy distribution (e.g. Doroshkevich 
et al. 1999). The most stringent constraints on cosmological 
models should be provided by direct observation of this evo- 
lution through measuring the properties of large-scale struc- 
tures at the highest possible redshifts; a task beyond the scope 
of the current large galaxy redshift surveys which reach only 
z ~ 0.2-0.3. One approach to investigate large-scale structure 
at higher redshifts is to use quasars, which are much eas- 
ier to detect at z > 1 than galaxies. At low-redshifts (z < 0.3), 
groups of quasars and AGN are shown to delineate the under- 
lying large-scale structure of galaxies and galaxy clusters (e.g. 
Sochting et al. 2002), a particular example being the "Great 
Wall" which is traced by a group of 19 AGN (Longo 1991). 
At higher redshifts (z > 0.5) a number of quasar superstructures 
have been discovered, made up of 10-25 quasars, and span- 
ning 50-200 /T ! Mpc (e.g. Webster 1982; Crampton et al. 1987, 
1989; Clowes & Campusano 1991, 1994). These large quasar 
groups (LQGs), if proved to be marking out volumes with a 
generalised galaxy overdensity, may represent the precursors 
of the "Great- Wall"-like structures such as those found in the 
local universe (Komberg & Lukash 1994), and hence may pro- 
vide ideal laboratories for studying the evolution of super-large 
scale structure of galaxies. 

The relationship between the quasar distribution and the 
underlying galaxy or mass distribution is likely to be both 
complex and redshift-dependent, particularly on cluster length- 
scales (<2/z _1 Mpc) where the dependence on the environ- 
ment of the fuelling mechanism of the quasar is likely to 
be dominant. On larger scales, the quasar spatial distribu- 



tion should mirror that of their host galaxies, the majority of 
which, from high-resolution observations of bright quasars, 
appear to be massive ellipticals (e.g. McLure et al. 1999; 
Kukula et al. 2001), which are known to be already in place 
at z ~ 1.2 (Blakeslee et al. 2003). This result is consistent 
with the requirement of a massive black-hole to power the 
quasar, in conjunction with the observed correlation between 
the masses of the black-hole and the bulge of the host galaxy 
(Merritt & Ferrarese 2001). The expectation that quasars and 
their host galaxies (i.e. massive ellipticals) trace the underly- 
ing mass distribution in similar ways also appears confirmed 
through comparison of the quasar and galaxy power-spectra 
(Hoyle et al. 2002) and 2-point correlation-functions (Croom 
et al. 2001) which show the same general form over the range 
2-30/T 1 Mpc. 

Of the quasar groups identified to date, the largest spa- 
tially is that of the Clowes-Campusano LQG, with at least 18 
quasars at 1.2 < z < 1.4 towards ESO/SERC field 927 (Clowes 
& Campusano 1991, 1994; Graham et al. 1995; Newman et al. 
1998; Clowes et al. 1999), with a spatial overdensity of 6-10 
times over the mean. The structure spans ~ 2.5° x 5° on the 
sky, corresponding to ~ 120 x 240/T 2 Mpc 2 atz ^ 1.3, making 
it the largest known structure in the high redshift universe. 

To investigate the underlying galaxy distribution traced by 
the Clowes-Campusano LQG, a study of Mg n absorption sys- 
tems in the spectra of background quasars has been carried out 
(Williger et al. 2002). Spectra for 23 quasars (1.23 < z < 2.68) 
in a 2.5° x 2.5° field towards the centre of the LQG were ob- 
tained. The advantage of using Mg n absorbers is that they trace 
much lower mass overdensities than quasars, and thus offer a 
much more detailed picture of the overall mass distribution. 
38 Mg ii absorbers were identified (Wo > 0.3A) while only 24 
would be expected given the quasars observed and the wave- 
length coverage of the survey. Of these, 1 1 were found to be 
associated with the LQG having 1 .2 < z < 1 .4 whereas only 
3.8 ±1.8 would be expected, a result significant at the 1-2% 
level and confirming the hypothesis that the quasar group delin- 
eates an underlying galaxy large-scale structure. The Mgn ab- 
sorber survey also identified a second structure at z ^ 0.8 made 
up of 7 Mg ii absorption systems at 0.77 < z < 0.89, again sig- 
nificant at the 1-2% level. This structure coincides with a group 
of 14 quasars at z ~ 0.8 previously identified by the Chile-UK 
Quasar Survey (CUQS, Newman 1999) as a structure signifi- 
cant at the 6.6% level. This structure spans ~ 3.5° x 3° on the 
sky, corresponding to ~ 120 x 100 /T 2 Mpc 2 at z ^ 0.8. 

Figure 1 shows the Clowes-Campusano LQG 
(left; 1 .2 < z < 1 .4) and the z =* 0.8 structure (right; 
0.77 < z < 0.93). Each quasar and Mg ii absorber within 
the redshift range of the two structures are indicated by 
stars and circles respectively, labelled by their redshift. The 
boundaries of the plots match the boundaries of the AQD 
survey of ESO/SERC field 927 (Clowes & Campusano 1991, 
1994), while the boundaries of the Mg ii survey (Williger et al. 
2002) and the BTC field (a 40 x 35 arcmin 2 region covered by 
deep V and / imaging; Haines 2001) are indicated by dashed 
and solid boxes respectively. The detected Mg n absorbers at 
z ^ 0.8 and z- 1.2 are found in the same regions as the large 
quasar concentrations detected at the same redshifts. 
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Fig. 1. The Clowes-Campusano LQG (left; 1.2 < z < 1.4) and the z - 0.8 structure (right; 0.77 < z < 0.93). Each quasar and Mgn absorber 
within the redshift range of the two structures are indicated by stars and circles respectively, labelled by their redshift. Double circles indicate 
two Mg ii systems towards a line of sight. The boundaries of the plots match the boundaries of the AQD survey of ESO/SERC field 927 (Clowes 
& Campusano 1991, 1994), while the boundaries of the Mgn survey (Williger et al. 2002) and the BTC field (Haines 2001) are indicated by 
dashed and solid boxes. 



To directly detect part of the underlying galaxy distribu- 
tion accompanying the large-scale structures observed in the 
quasar and Mg n absorber distributions at z - 0.8 and z - 1.2, 
we have also undertaken an ultra-deep optical imaging sur- 
vey of a 40 x 35 arcmin 2 region - the BTC field - containing 
three quasars from the Clowes-Campusano LQG, in both V 
and / passbands. These data were used previously by Haines 
et al. (2001b) in conjunction with UKIRT /f-band imaging 
to identify a cluster at z - 1.2 associated with the z = 1.233 
LQG quasar (redshift value changed from z = 1.226 as re- 
ported in Haines et al. 2001b based on new spectroscopic data; 
Williger et al. 2002). A factor -11 overdensity of / - K > 3.75 
galaxies was identified in a 2.25 x 2.25 arcmin 2 field centred 
on the z = 1.233 LQG quasar, for which NIR imaging had 
been obtained, reaching K ^ 20. In particular, 15-18 galaxies 
with colours consistent with being a population of passively- 
evolving massive ellipticals at the quasar redshift were found. 
This field was reobserved using the New Technology Telescope 
(NTT) with its infrared camera (SOFI) for the study described 
in this article, and is part of field 1 (see Sect. 5). At the 
redshift of the Clowes-Campusano LQG (z = 1.2) the BTC 
field covers 31 x 27 /T 2 Mpc 2 , while at z = 0.8 the field cov- 
ers 23 x 20/T 2 Mpc 2 . These observations along with follow-up 
NIR imaging are described in Sect. 2. 

To identify any possible clustering or large-scale structure 
in the form of galaxies at z ~ la variant of the cluster red 
sequence (CRS) algorithm of Gladders & Yee (2000) is ap- 
plied to the BTC dataset. This approach is motivated by the 
observation that the bulk of early-type galaxies in clusters lie 
along a tight, linear colour-magnitude relation, the CRS, that 
becomes increasingly red with redshift (depending on choice 
of filters). In Sect. 3 we describe the variant of the CRS al- 
gorithm used, whereby galaxies redder than the z = 0.5 CRS 



are selected, and demonstrate its ability to identify early-type 
galaxies at 0.5 <z< 1.3, i.e. the redshift range of interest. In 
Sect. 4 the number density of these red galaxies in the BTC 
field is compared with that of a similar control field (in terms 
of passbands, depth and area covered) taken from the Deep 
Lens Survey (DLS) of Wittman et al. (2002). Evidence is found 
for there being one or more large-scale structures at z ~ 1 in 
the BTC field manifested as an overall excess of red galaxies. 
In Sect. 5 we describe follow-up targeted deep near-infrared 
imaging of four z ~ 1 cluster candidates identified as density 
enhancements in the red galaxy spatial distribution which also 
have apparent red sequences at V - I - 2.8. The implications 
of our results are discussed in Sect. 6 and the conclusions pre- 
sented in Sect. 7. 

Throughout the paper we assume a flat, A-dominated cos- 
mology with Q. m = 0.3, Q.a = 0.7, and for the predictions of the 
evolution of galaxy colours we assume H„ = 70kms" 1 Mpc" 1 . 
With these parameters, the age of the Universe is 13.5 Gyr, and 
the redshifts z = 0.8 and z — 1 .2 correspond to look-back times 
of 6.8 and 8.4 Gyr respectively. 



2. OBSERVATIONS 

2.1. Optical imaging data 

The optical data were obtained using the Big Throughput 
Camera (BTC) on the 4-m Blanco telescope at the Cerro 
Tololo Inter- American Observatory on April 21/22 and 22/23 
1998. The resulting BTC field (as we shall refer to it) covers 
40 x 35 arcmin 2 in both V and / imaging, reaching V - 26.35 
and / * 25.85. It is centred at 10H7 ffl 30 5 , +05°32'00" (J2000) 
and was selected to contain three members of the Clowes- 
Campusano LQG. 
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The BTC is made up of four 2048 x 2048 CCDs which 
have pixels of size 0.43 arcsec, giving a field of view for each 
CCD of 14.7 x 14.7 arcmin 2 . The CCDs are arranged in a 2 by 
2 grid and are separated by gaps of 5.6 arcmin. To obtain a con- 
tiguous image the camera had to be shifted between exposures. 
The camera geometry meant that it was not possible to obtain 
a uniform coverage of the field. This, along with one of the 
CCDs being markedly less sensitive (55% less in V, 30% in 
/) than the other three, meant that there were significant varia- 
tions in depth across the field, and so separate magnitude limits 
have been determined for each subregion discussed in this pa- 
per. Further details of the observing and reduction procedures 
for the BTC images are presented elsewhere (Haines 2001). 



2.2. Near-infrared imaging data 

Deep near-infrared imaging data were obtained using the SOFI 
camera on the 3.5-m ESO New Technology Telescope (NTT) 
on April 1-2 2002. The SOFI camera has a Rockwell HgCdTe 
1024 x 1024 Hawaii array with a pixel size of 0.292 arcsec 
giving a field of view of 4.9 x 4.9 arcmin 2 which corresponds 
to comoving angular scales of 3.8 x 3.8 h~ 2 Mpc 2 at z = 1.2 
and 2.8 x 2.8 /r 2 Mpc 2 at z = 0.8. Given the field of view, tar- 
geted fields within the BTC field were selected as having 
high-densities of red galaxies and apparent red sequences at 
V - I - 2.8 that may indicate galaxy clusters at z ~ 1. The best 
three candidates were observed for 3600s through both J and 
K s filters reaching J ^ 23 and K s ^ 21, and another four fields 
were observed for 900s in K s reaching K s ^ 20. It was found 
after the first night that the sky background levels in the K s - 
band fell by a factor of three during the first half of the night, 
and so field 1, which had been observed at the start of night one 
in K s , was reobserved for a further 3600s on the second night 
in more favourable conditions. Conditions were photometric 
throughout and the seeing was 0.9-1.5 arcsec. 

The NIR imaging data were reduced using standard IRAF 
routines. To allow the estimation of sky levels to be made from 
the observations themselves, each integration of 60s duration, 
and made up of 3(6) sub-integrations for the J(K s )-bands, was 
jittered randomly in a 20 arcsec wide box. For each individ- 
ual exposure, the sky background was estimated by combin- 
ing the 10 integrations closest in time by the median, after re- 
jecting all pixel values greater than 3 <x from the median. The 
sky image is then subtracted from its corresponding image af- 
ter being scaled to have the same median value. The images 
are then flat-fielded through the use of dome flats taken using 
the SOFIJma_calJ>pecialDomeFlat template which allows a 
residual shade pattern to be estimated and removed. An illu- 
mination correction was then applied to the images using illu- 
mination correction surfaces produced by observing a photo- 
metric standard star at differing parts of the array and fitting a 
surface. The SOFI array suffers from interquadrant row cross 
talk, which produces ghosts that affect all the rows of a bright 
source, and also the corresponding rows in the opposing half of 
the detector. Although the effect is not completely understood, 
it is well described and was corrected for in each of the raw 



exposures. The individual exposures were then registered and 
coadded to produce the final images. 

3. SELECTION OF EARLY-TYPE GALAXIES AT 

0.5 < z < 1.3 

To detect and identify any possible large-scale structure at z ~ 1 
a variant of the cluster red sequence algorithm of Gladders & 
Yee (2000) is applied to the BTC dataset. This approach is mo- 
tivated by the observation that the bulk of early-type galaxies 
in all rich clusters lie along a tight, linear colour-magnitude (C- 
M) relation - the cluster red sequence. Studies of low-redshift 
clusters (e.g. Bower et al. 1992) indicate that, irrespective of 
the richness or morphology of the cluster, all clusters have red 
sequences, whose k-corrected slopes, scatters and colours are 
indistinguishable. This indicates that the early-type galaxies 
which make up the red sequences form a homogeneous pop- 
ulation, not only within each cluster, but from cluster to clus- 
ter, and also that red sequences are universal and homogeneous 
features of galaxy clusters, at least at z < 0.2. The photomet- 
ric evolution of the cluster red sequence with redshift has been 
studied (e.g. Aragon-Salamanca et al. 1993; Ellis et al. 1997; 
Stanford et al. 1998; Kodama et al. 1998) for clusters out to 
z - 1.2 indicating: that red sequence remains a universal fea- 
ture of clusters; that the stellar populations of its constituent 
early-type galaxies are formed in a single, short burst at an early 
epoch (z/>3); and that the galaxies have evolved passively 
ever since. This is confirmed by the spectra of red-sequence 
galaxies which in nearby clusters are best fit by simple stellar 
populations of ages 9-12 Gyr, resulting in their having charac- 
teristically strong spectral breaks at 4000 A, and correspond- 
ingly red U - V colours. 

The basic implementation of the cluster red sequence algo- 
rithm is to split the galaxy catalogue into a series of redshift 
slices by selecting the subset of galaxies whose colours and 
magnitudes are consistent with the C-M relation of clusters at 
that redshift. For each slice the galaxy surface density is then 
estimated, and significant peaks in each slice are then identi- 
fied as cluster candidates at that redshift. The key advantages 
of this algorithm for cluster detection are: that it requires only 
photometry in two passbands, as is the case here; the colour 
of the observed red sequence can be used as a precise red- 
shift indicator; and it does not suffer from projection effects 
or foreground contamination. The latter two are contingent on 
the passbands used covering the rest-frame 4000A break for 
the redshift range of interest, as then the evolution of colour 
of the red sequence with redshift is greatest allowing the most 
accurate results, and also the red sequence galaxies are as red 
or redder than all galaxies at that redshift and all nearer galax- 
ies, hence eradicating the problem of foreground contamina- 
tion. Hence for our data with V and / photometry, the optimal 
redshift range is 0.3 <z < 0.8, and conversely the optimal pass- 
bands for estimating redshifts of clusters at 0.8 <z< 1.3 are / 
and K s . 

Figure 2 shows the evolution of galaxies' V-I colours as a 
function of galaxy type and redshift. The solid curves indicate 
the expected evolution of the red sequence with redshift, the 
most massive red sequence galaxies being more likely to form 
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Fig. 2. Evolution of galaxies' V-I colours as a function of galaxy type 
and redshift. The solid curves represent stellar populations formed in 
an instantaneous burst at z/ = 3.5 (thick) and Zf = 2.0 (thin), and 
are thought to represent the colour evolution of massive ellipticals. 
The dashed and dot-dashed curves represent stellar populations with 
exponentially-decaying star-formation rates with time-scales (t) of 1 
and 5 Gyr respectively, and are thought to represent the colour evolu- 
tion of disk-dominated galaxies. The points correspond to the colour 
of the red sequence for clusters from the surveys of SED98, S97, A93 
and N01. The dotted line indicates the colour selection threshold for 
red galaxies at / = 23.5. 



early (z/ - 3.5) as shown by the thick line, while the less mas- 
sive early-type galaxies are more likely to follow the thin line, 
indicating later formation epochs (zf - 2.0) or lower metallic - 
ities. The dashed and dot-dashed curves correspond to stellar 
populations with exponentially-decaying star-formation rates 
with time-scales (t) of 1 and 5 Gyr respectively, and are thought 
to represent the colour evolution of disk-dominated galaxies. 
The points correspond to the V - I colour of the red sequence 
for clusters from the surveys of Stanford et al. 1998 (SED98 : 
6 clusters at 0.54 < z < 0.7), Smail et al. 1997 (S97 : 3 clusters 
at z ~ 0.54), Aragon-Salamanca et al. 1993 (A93 : 10 clusters 
at 0.5 < z < 0.9) and Nelson et al. 2001 (N01 : 30 clusters at 
0.3 < z < 0.9). The figure indicates that at redshifts less than 
Z ~ 0.8 the V - I colour of the red sequence increases monoton- 
ically with redshift, and so can be used to efficiently estimate 
the cluster redshift. At higher redshifts the V-I colour becomes 
a less efficient redshift indicator, as its dependence on redshift 
decreases, and dependence on star-formation history increases. 
The V-I colour of the red sequence in clusters at 0.8 < z < 1 .3 
remains roughly constant at V - I - 2.8, as the 4000A break 
passes through the /-band. However the two models used to 
describe the red sequence galaxies show significant divergence 
at z > 1, indicating that the colour of red sequences in clus- 
ters at 0.8 < z < 1.3 could be in the range 2.5 < V - 1 < 3.2, 
and so it is not possible to obtain an accurate redshift esti- 
mate based on the V-I colour of the red sequence for z ~ 1 
clusters. At these redshifts the V-band now corresponds to a 



rest-frame wavelength of 2500-3000A, so the V - I colour is 
affected by even small amounts of recent star-formation as is 
indicated by the divergence of the exponentially -decaying star- 
formation rate models from the burst models. A third problem 
with the basic cluster red sequence algorithm at these redshifts 
is that at z > 1 even L* galaxies are reaching the detection limit 
in V, and so that the red sequence is likely to be smeared in the 
colour direction by the less accurate photometry, reducing the 
signal of any cluster. 

Given all these problems facing the standard cluster red se- 
quence algorithm, it was decided that instead of considering 
narrow redshift slices over the redshift range 0.8 <z< 1.3, a 
single high-redshift slice would be used to identify structures 
at this redshift range, containing all galaxies redder in V - I 
than the cluster red sequence at z — 0.5. Kodama et al. (1998) 
show that the Kodama & Arimoto (1997) evolutionary model 
for elliptical galaxies with zj =4.5 describes well the evolu- 
tion of the zero-point and slope of the cluster red sequence over 
0.3 < z < 1 -2. We take directly from Figs. 4 and 5 of Kodama et 
al. (1998) their values for the zero-point and slope of the clus- 
ter red sequence at z = 0.5, resulting in the selection criterion 
for red galaxies of 



V-I < 3.689 - 0.0652 x W; 



'total 



<23.5. 



(1) 



A magnitude limit is also applied, as most red sequence galax- 
ies should be brighter than / = 23.5 to z - 1.3 (an L* early- 
type galaxy has / ^ 23.0 at z = 1.3), and any galaxy with a 
sufficiently red V-I colour will be close to the V-band de- 
tection limit. Hence all galaxies classed as belonging to the 
high-redshift slice will have V - I > 2.15, as indicated in Fig. 2 
by the horizontal dotted line. An examination of this fig- 
ure indicates that this approach will remove all foreground 
galaxies (z < 0.5), as well as any star-forming galaxies, leav- 
ing only passively-evolving galaxies in the redshift range of 
0-5 <z< 1.3. It should then be possible to identify any galaxy 
clusters or large-scale structures at this redshift range as den- 
sity enhancements in the spatial distribution of these red galax- 
ies (as we shall refer to them), but there will be little ability to 
distinguish across this redshift range. 

3.1. Red Galaxies in the ESO Imaging Survey of the 
Hubble Deep Field South 

As a test of the efficiency of the colour-magnitude criteria to 
select early-type galaxies at 0.5 <z< 1.3, we have applied the 
same criteria to galaxies in the ESO Imaging Survey (EIS) of 
the Hubble Deep Field South. This comprises deep optical data 
for a 25 arcmin 2 field reaching 2<x limiting magnitudes of U ~ 
26, B ~ 26.5, V ~ 26, R ~ 26, / ~ 24.5 obtained using the 
SUSI2 imager on the NTT, plus near-infrared data reaching / ~ 
24, H ~ 22.5 and K s ~ 22 for about 40% of the optical field. 
In total 34 red galaxies (those satisfying Eq. 1) are identified in 
the EIS-DEEP field, of which 9 have near-infrared imaging. 

To estimate the photometric redshifts of each of the galax- 
ies from its UBVRI (and where possible JHK S ) photometry, 
and to produce the model colours tracks used throughout this 
paper, the hyperz code of Bolzonella et al. (2000) has been 
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Fig. 3. Best-fit synthetic SED to sample galaxy in EIS-DEEP field 
with UBVRIJK, photometry 



used. It builds synthetic template galaxies using the Bruzual 
& Chariot evolutionary code (GISSEL98; Bruzual & Chariot 
1993). It has stellar populations with eight star-formation his- 
tories roughly matching the observed properties of local galax- 
ies from E to Im type: an instantaneous burst; a constant star- 
formation rate; and six exponentially-decaying star-formation 
rates with time-scales, t, from 1 to 30 Gyr chosen to match the 
sequence of colours from E-S0 to Sd galaxies. The models as- 
sume solar metallicity and a Miller & Scalo (1979) initial mass 
function with upper and lower mass-limits for star-formation 
of 125M and 0.1 M G respectively. Internal reddening is also 
considered through the Calzetti et al. (2000) model with A v 
allowed to vary between and 1 mag. The hyperz software 
then produces a photometric redshift probability distribution 
for each galaxy through a chi-squared minimization process, 
allowing for all galaxy ages, star-formation histories and Ays. 
Hence instead of a single best-fitting redshift for a galaxy, a 
range of compatible redshifts is produced. 

We find that, despite allowing for all galaxy types and ages, 
all the red galaxies are best-fit by early-type galaxies whose 
stellar populations formed in an instantaneous burst. As pre- 
dicted, the selection of galaxies redder in V - 1 than the z = 0.5 
cluster red sequence removed all foreground contamination, 32 
of the 34 red galaxies in the EIS-DEEP field are constrained 
at the lcr level to be at z > 0.6, the other two being con- 
strained to z > 0.45. The selection of red galaxies effectively 
identifies a population of passively-evolving early-type galax- 
ies at 0.5 <z< 1.3, whose strong 4000A breaks are responsi- 
ble for their red V - I colour. This is demonstrated by Fig. 3 
which shows the best-fitting synthetic SED to a sample galaxy 
(or=22:32:42.31, £=-60:33:43.0) in the EIS-DEEP field with 
UBVRIJK S photometry. The redshift estimate is constrained 
to a narrow range with Zbest = 0-962+^^, and it is clear that 



this is due to the very strong 4000A break which is manifested 
as the red V - I colour. 

4. RESULTS : THE BTC FIELD 

4.1. Red galaxies in the BTC field 

To identify the red galaxies in the BTC field, catalogues of 
sources in the registered V and / images were created, using 
the / image to detect the sources given the red colour of the 
sources of interest. Registration and astrometric calibration of 
the images was performed using the USNO astrometric cata- 
logue, and is accurate to ^ 0.05 arcsec. Photometric calibration 
of the V and / images onto the Johnson-Kron-Cousins system 
was obtained using 15 Landolt (1992) standard stars observed 
at varying airmasses, resulting in zero-point uncertainties of 
0.026 (V) and 0.009 (/), producing an uncertainty in V - 1 of 
0.028 mag. The V-band uncertainty is greater due to the large 
differences in sensitivity across the four CCDs. Object detec- 
tion was performed on the /-band image using SExtractor 
(Bertin & Arnouts 1996) for objects with 4 contiguous pix- 
els 0.8cr over the background level. The total / magnitude was 
taken to be the MAG_BEST output from SExtractor and the 
V - I colour determined using fixed apertures of 2.5 arcsec di- 
ameter with SExtractor in two image mode. 

All sources satisfying the red galaxy selection criteria (Eq. 
1) were then visually verified in the / image, and spurious 
objects removed. The bulk of these contaminant sources are 
due to random noise events near the field edges, and given the 
faint magnitude limit of / = 23.5 is 2 mags above the detection 
threshold of the /-band images, it is relatively easy to discrim- 
inate between real and spurious sources in the vast majority of 
cases. Regions around bright, heavily-saturated stars were also 
rejected as in the / image the wide wings of the point-spread- 
functions artificially boost the /-band flux of faint sources, re- 
sulting in rings of "red" sources. In total ~ 5% of the area of 
the BTC field is discarded in this manner, but as stars should 
be randomly distributed this should not affect the results. 

The most significant source of contamination of the red 
galaxy subset is that from faint red dwarfs. To minimise this 
contamination, the in-built star-galaxy classifier of SExtractor 
was applied, and 467 sources with CLASS_STAR > 0.95 taken 
to be stars and rejected. The stellarity classifier is believed to be 
efficient to / ~ 22 where the vast majority (> 90%) of sources 
have stellarities greater than 0.95 (i.e. stars) or less than 0.2 
and are hence galaxies. At fainter magnitude levels there is 
likely to remain some contamination due to stellar sources, as 
the classifier becomes less efficient, and although at this level 
there are two clear loci in the CLASS-STAR distribution, pre- 
sumably corresponding to sources well-classified as stars or 
galaxies, there remain ~ 250 sources with ambiguous stellar- 
ities (0.2 < CLASS. STAR < 0.95). 

After removing spurious and stellar sources, and discarding 
the regions whose photometry is affected by bright stars, 2674 
red galaxies remain in a total area of 1277 arcmin 2 . This gives a 
density of 2.09 ± 0. 14 galaxies arcmin -2 , where the uncertainty 
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considers Poisson noise as well as the effect of the two-point 
angular correlation function co(0) through the equation 



a 2 



n + n 



_ 2 J / oj(0 12 )da l da 2 

/ JdQ.idQ.2 



(2) 



where h is the number of galaxies, and a>{6\2) = A" d (l') -1 - 045 
is the two-point angular correlation function of red galaxies 
in the BTC field (Haines 2001). The red galaxies are much 
more strongly clustered than galaxies selected by magnitude 
only (i.e. / < 23.5), with A^ d (l') = 0.0863 in comparison to 
A^ <235 (T) = 0.0175. As a consequence, the effect of the clus- 
tering on the rms uncertainty in galaxy counts is significant, 
three times that of Poisson noise alone. 

4.2. An excess of red galaxies - A large-scale 
structure atz ~ 1 across the BTC field 

Any large-scale structure at z ~ 1 underlying either the Clowes- 
Campusano LQG at z - 1 .2 or the structure of quasars and 
Mgn absorbers at z - 0.8 should manifest itself as an excess 
of red galaxies across the whole BTC field. To determine the 
extent of any excess we estimate the expected field density of 
red galaxies through consideration of a suitable random control 
field. For this purpose we have used field Flp22 from the Deep 
Lens Survey (DLS) of Wittman et al (2002; http://dls.bell- 
labs.com). The Deep Lens Survey is an ongoing ultra-deep 
multi-band optical survey planned to cover seven randomly- 
chosen 4deg 2 fields with deep BVRz imaging obtained us- 
ing the Mosaic cameras on the 4-m telescopes at the Cerro 
Tololo and Kitt Peak observatories. The publicly-available field 
Flp22 has BVRIz imaging covering 35 x 35 arcmin 2 centred 
at h 53 m 25 s .30, +12°33'55'.'0 and has limiting magnitudes of 
/ ~ 25 and V ~ 26. It is an excellent comparison field, covering 
a similar size, and reaching similar depths in the same filters as 
the BTC field. 

To identify the red galaxies in the DLS data, exactly the 
same detection and source extraction processes are applied as 
for the BTC field, with objects detected in the / image using 
SExtractor and colours determined through 2.5 arcsec diam- 
eter apertures using SExtractor in two-image mode. As for 
the BTC field, red galaxies are identified as sources redder in 
V-I than the z = 0.5 cluster red sequence and brighter than 
/ = 23.5. Once again, regions around bright stars are dis- 
carded, spurious objects are rejected through visual inspec- 
tion, and stars identified by the star-galaxy classifier as sources 
with CLASS_STAR > 0.95 and removed from the final cata- 
logue. We expect the efficiency of the identification and re- 
moval of both stars and spurious objects from the red galaxy 
catalogues to be similar for the BTC and DLS data, as both 
have comparable magnitude limits (/ « 25.5) and levels of see- 
ing (« 1.2 arcsec), resulting in near-identical image qualities. 

In total 1953 red galaxies are found in the 1181 arcmin 2 
area of the DLS image not affected by bright stars, correspond- 
ing to a mean density of 1.65 ± 0.12 gals arcmin -2 . 

We would thus expect 2111 ± 151 galaxies in the BTC 
field, where the uncertainty assumes both Poisson noise and the 
effect of the two-point angular correlation function as before. 




V-I 



Fig. 4. The V-I colour distributions of the observed excess of red 
galaxies for the BTC field (solid curve connecting filled circles). For 
comparison the dotted curve shows the colour distribution of galax- 
ies in the DLS field (scaled to fit plot), and the star-shaped symbols 
indicate the colour distribution of stars in the BTC field. 



The density of red galaxies in the BTC field is 27% greater than 
that observed for the DLS image, corresponding to an excess of 
563 galaxies over the 40 x 35 arcmin 2 field, a result significant 
at the 3.7 cr level. 

To ascertain the nature of this excess we plot in Fig. 4 
the V - I colour distribution of the galaxy excess in 0.1 mag 
wide bins in V - I as filled circles connected by a solid 
line. This is determined by scaling the colour distribution of 
galaxies in the DLS to have the same effective area as the 
BTC field, and subtracting it from the colour distribution of 
galaxies in the BTC field. The errors shown are determined 
as cr 2 = 7V,(BTC) + Af,(DLS). Significant excesses are appar- 
ent for all bins in the range 2.5 < V - 1 < 3.8. The excess 
shows a coherent peak at 2.7 < V - I < 3.1, and is greatest 
for the 2.8 < V - 1 < 2.9 bin for which 224 such galaxies are 
identified in the BTC field whereas only 135 would be ex- 
pected given the density in the DLS image, corresponding to 
a 66% excess, significant at the 3.5 cr level. A comparison 
with Fig. 2 shows that this peak corresponds to the colours ex- 
pected of passively-evolving early-type galaxies at the redshift 
range of interest, 0.8 <z< 1.3. For comparison the star sym- 
bols indicate the colour distribution of stars (defined to have 
CLASS_STAR > 0.95) in the BTC field. The colour distribu- 
tions of stars and the excess red galaxies appear quite different, 
indicating that stellar contamination is not a major contributor 
to the red galaxy excess. The colour distribution of galaxies in 
the DLS field (scaled to fit plot) is also shown by the dotted 
curve, and again appears different to the distribution of the ex- 
cess red galaxies, indicating that the excess forms a population 
that is quite distinct to that from the field. 

Figure 5 shows the excess of red galaxies binned by / mag- 
nitude as a solid curve connecting filled circles. This is deter- 
mined by scaling the magnitude distribution of galaxies in the 
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19 20 21 22 23 

I 

Fig. 5. The / magnitude distribution of the observed excess of red 
galaxies across the BTC field shown as a solid curve connecting filled 
circles. The dashed curve corresponds to the best-fitting Schechter 
function with a = -0.3 and m"j = 20.90. The dotted curve shows 
the magnitude distribution of the 563 galaxies nearest the z = 0.5 
cluster red sequence in the BTC data, while the star symbols indicate 
the magnitude distribution of stars redder than the z = 0.5 cluster red 
sequence. 

DLS image to have the same effective area as the BTC field, 
and subtracting it from the magnitude distribution of red galax- 
ies in the BTC field. Significant excesses are observed for all 
bins in the range 20 < l tot < 23.5. The excess is particularly 
significant for 20 < I, , < 21 where 364 red galaxies are found 
in the BTC field whereas only 209 would be expected from the 
DLS data, corresponding to an excess of 74%, and a signifi- 
cance level of 4.3 <x, assuming the same effect on statistics of 
galaxy clustering as previously. 

For comparison the magnitude distribution of stars redder 
than the z = 0.5 cluster red sequence (i.e. having the same 
colours as the red galaxies) in the BTC field is indicated by the 
star symbols. The magnitude distribution of red stars appears 
flat, and is quite distinct from that of the red galaxy excess, 
again confirming that the red galaxy excess is not primarily 
due to stellar contamination. In particular, at bright magnitudes 
(7 < 21) where the stellar classifier is virtually 100% efficient, 
the two magnitude distributions are completely different. 

To examine whether the excess could be due to a large 
photometric error producing a colour shift in either the BTC 
or DLS data, we introduce an artificial shift in the V - 1 
colours of the BTC galaxies sufficiently large to remove the 
observed excess in the BTC field. To account for the whole ex- 
cess in this manner requires a V - I colour shift of -0.15 mag. 
Assuming zeropoint calibration uncertainties of 0.01 mag for 
both V and / DLS images, we obtain an overall relative pho- 
tometric calibration uncertainty between the BTC and DLS 
data of A(V - I) = 0.03. Hence, a calibration error as large as 
0.15 mag appears unlikely, but as a further check we consider 
the magnitude distribution of those galaxies which are initially 



classed as red galaxies, but which would be lost if there was a 
V — I colour shift of -0.15 mag, i.e. those with 

3.689 - 0.0652 x / < V - I < 3.839 - 0.0652 x / ; 
Uotal < 23.5. 

The magnitude distribution of these galaxies is shown in Fig. 5 
as a dotted curve, and is clearly substantially different to the 
magnitude distribution of the observed excess, the former be- 
ing unable to reproduce the sharp increase at / ^ 20 observed 
in the latter. The distinctness of both the colour and magnitude 
distributions of the observed red galaxy excess in comparison 
to those predicted for excesses produced by stellar contami- 
nation or photometric calibration errors, in particular the peak 
at V - 1 2.8 and the sharp increase at / 20, indicates that 
the red galaxy excess observed is not primarily due to either of 
these sources. 

Given the relative ease of differentiating between real and 
spurious sources in the image, that although their removal is a 
subjective process, we estimate the uncertainty resulting from 
this process to be at the 1% level, or ~ 25. The major uncer- 
tainties in the final figure are the level of stellar contamina- 
tion, which could contribute as many as ~ 250 sources, and 
the uncertainty in the zero-point of the V - I colour, which we 
take to be A(V - I) = 0.03. resulting in an uncertainty in the 
red galaxy numbers of 1 17 (the number of red galaxies which 
would be lost if their V - I colour were reduced by 0.03 mag.) 
The effect of uncertainty in the /-band magnitude is much less, 
even assuming a pessimistic level of Al = 0.02 produces an un- 
certainty of only 18 red galaxies (the number of galaxies with 
23.48 < / < 23.50). 

The galaxy luminosity function is a powerful tool in cos- 
mology by allowing the bulk properties of galaxies to be mea- 
sured. Numerous studies show that the luminosity function is 
well described by the Schechter (1976) function 

<KL)dL = f(±) a e-^d(±), (4) 
which is written as a function of apparent magnitude as 

<f>(m)dm = kfe [ - kla+l)(m - m ' ) - e ~ k{ """' ) 1 dm, (5) 

where k= 0.4 In 10. 

The redshift evolution of the galaxy luminosity function 
has been examined for cluster populations out to z ~ 1 (e.g. de 
Propris et al. 1999), finding that L* evolves in a manner con- 
sistent with the passive evolution of galaxies formed at early 
epochs, although there remains some debate over the effect of 
galaxy mergers on <p{L). Studies measuring the evolution of 
m* with redshift for galaxy clusters to z ~ 1 (e.g. Nelson et 
al. 2001; de Propris et al. 1999) indicate that m*(z) can be fit 
by a single function with intrinsic scatter between clusters at 
the same redshift of cr(m*) ~ 0.2. This then affords the pos- 
sibility of estimating the redshift of a cluster from its galaxy 
luminosity function, and Nelson et al. (2001) indicate that 
the rms error of redshift estimates from m*(z) for 44 clusters 
at 0.3 < z < 0.9 is A(z) = 0.06. By calculating the luminosity 
function of ~ 75 000 galaxies from the 2dF Galaxy Redshift 
Survey for different subsets defined by their spectral proper- 
ties, Madgwick et al. (2002) show that a varies significantly 
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Fig. 6. Contour plot of the estimated density distribution of red galaxies in the BTC field (corresponding to 31 x 27 /T 2 Mpc 2 at z = 1.2). The 
dashed contour corresponds to the density of 1.65 galaxies arcmkr 2 , that observed for the DLS data. The shaded regions indicate the density 
peaks with solid contours corresponding to overdensities of 100%, 200% and 300% in comparison to the DLS data. The four quasars at z ~ 1.3 
are indicated by star symbols labelled by redshift. The NTT fields that targeted the four best cluster candidates at z ~ 1 are indicated by the 
labelled boxes. 



with spectral type. They find a systematic steepening of the 
faint-end slope, a, from passive (a = -0.54 ± 0.02) to active 
star-forming (a = -1.50 ± 0.03) galaxies. 

We fit the Schechter function to the observed excess 
through a chi-squared minimization process, allowing m* n (p* 
and a to vary. The best-fitting function is shown as the dashed- 
curve in Fig. 5 and the optimal values for the free parameters 
found to be m] = 20.90 ± 0.12 and a = -0.3 ± 0.25. 

The observed shallow faint-end slope is comparable to that 
for local early-type galaxies (Madgwick et al. 2002), and hence 
confirms previous suggestions that the excess is made up of lu- 
minous early-type galaxies. It also appears to indicate that the 
excess is due to a structure at a single redshift, as a wide spread 
of redshifts would produce a steeper faint end curve as the lu- 
minosity function is convolved with the redshift distribution. 

In determining m*(z) Nelson et al. (2001) use a faint-end 
slope of a = - 1 .25, and to compare our data with theirs, and 



to estimate the redshift of the excess galaxies we find that for 
a fixed a = -1.25 we obtain m* = 20.12 ± 0.12, which results 
in a redshift estimate of z = 0.83 ± 0.08. It should be noted that 
the excess magnitude distribution we are measuring is colour 
selected and biased towards early-type galaxies, and this ex- 
plains the shallow faint-end slope obtained, whereas Nelson 
et al. (2001) consider the overall galaxy excess as a function 
of magnitude. However as early-type galaxies dominate the 
bright-end of the galaxy luminosity function, this should not 
affect the redshift estimate significantly. 

The fact that we are able to fit a Schechter function to the 
magnitude distribution of the excess, plus the observed peak 
in the colour distribution, indicates that this excess is predomi- 
nately due to a superstructure containing many luminous galax- 
ies at a single redshift of z - 0.8. This corresponds to the red- 
shift of one of the two structures identified in the quasar and 
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Fig. 7. The V - I against / colour-magnitude diagrams for galaxies in the four NTT fields. Galaxies in the circular regions corresponding 
to the clusters identified in each field are indicated by squares, with those having (not having) the colours expected of early-type galax- 
ies at the cluster redshift indicated by solid (open) symbols. For field 1 this corresponds to z — 1.2, and the galaxies are selected to have 
V - I > 1.8 and 3.6 < I - K s < 4.8, whereas for fields 2, 3 and 4 this corresponds to z - 0.8 and the galaxies are selected to have V - I > 2.3 
and 3.00 < I - K s < 3.75. The remaining galaxies in the fields are indicated by small crosses. 



Mgn absorber distributions, suggesting they are all manifesta- 
tions of the same underlying large-scale structure. 

4.3. Spatial distribution of the red galaxies - cluster 
candidates atz~ 1 

The overall excess of red galaxies across the BTC field is most 
likely due to the presence of one or more large-scale structures 
at z ~ 1 . As discussed previously, the most natural explanation 
is that either or both of the structures identified in the quasar 
and Mg n absorber distribution, the Clowes-Campusano LQG 
at z - 1.2 and the structure at z - 0.8, bely the presence of 
an underlying large-scale structure in the form of galaxies. To 
examine any relation between the red galaxy excess and the 
Clowes-Campusano LQG, and to identify suitable cluster can- 
didates for follow-up near-infrared studies, the density distri- 
bution of the red galaxies (i.e. those satisfying Eq. 1) across 
the BTC field is estimated through an adaptive kernel approach 
(Pisani 1993, 1996; Haines 2001), as shown in Fig. 6. The four 
quasars at z ~ 1.3 are indicated in the figure by star symbols 
labelled by redshift. 

Four good z ~ 1 cluster candidates are identified across 
the BTC field as the most significant density enhancements 
in the red galaxy spatial distribution, and also having appar- 



ent red sequences at V - 1 — 2.8, thus minimising the pos- 
sibility of projection effects. These were targetted for NIR 
imaging as shown in Fig. 6 by labelled boxes. Fig. 7 shows 
the V - 1 against / colour magnitude diagrams for galaxies 
in the NTT fields covering the four best z ~ 1 cluster candi- 
dates. Galaxies in the circular regions corresponding to the 
clusters identified in each field (see Fig. 9) are indicated by 
squares, with those having (not having) the colours expected of 
early-type galaxies at the cluster redshift (based on the VIK S 
colours) indicated by solid (open) symbols. For field 1 this 
corresponds to z — 1.2, and the galaxies are selected to have 

V - 1 > 1.8 and 3.6 < / - K s < 4.8, whereas for fields 2, 3 and 
4 this corresponds to z — 0.8 and the galaxies are selected to 
have V - I > 2.3 and 3.00 < I - K, < 3.75 (see Sect. 5.2). The 
remaining galaxies in the fields are indicated by small crosses. 
Each of the four NTT fields have excesses of 20-50 red galax- 
ies in comparison to the DLS data, and between them account 
for 140 excess red galaxies, or 25% of the total over the BTC 
field. Hence the clusters are not the sole reason for the excess, 
and the majority is due to a general large-scale excess. The 
only other cluster candidate (10 h 48 m 31 s , +05°23') with a com- 
parable excess of (^ 20) red galaxies has a red sequence at 

V — I ~ 2.4 suggesting it is at a lower redshift. 
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Table 1. Coordinates and magnitude limits for the fields covered by 
near-infrared imaging. 



5. NEAR INFRARED IMAGING 

Photometric calibration of the NIR data was performed through 
the repeated observations of the faint near-infrared NICMOS 
standard stars of Persson et al. (1998). In total 24 standards 
were observed through the 2 nights, matching the variation of 
airmass of the science exposures, resulting in zeropoints for 
each set of science data accurate to ~ 0.01 mag for both J and 
/f s -bands. The V-, /-, and /-band data were registered with 
the K s -bnnd data, and the images convolved with a Gaussian 
kernel to match the PSF of the image with the worst seeing. 
Object detection was performed on the 7T,-band image using 
SExtractor (Bertin & Arnouts 1996), and colours determined 
using fixed apertures of 2.5 arcsec diameter with SExtractor in 
two-image mode. Photometric uncertainties were determined 
from the noise levels in the pre-smoothed images. 

The coordinates and magnitude limits of each of the fields 
are presented in Table 1 . Field 1 corresponds to the galaxy clus- 
ter associated with the z = 1.233 LQG quasar, which had been 
previously identified from K imaging (Haines et al. 2001b). 
Fields 2 and 3 both contain double clusters identified as over- 
densities in the red galaxy spatial distribution in the BTC field. 
Field 4 contains a cluster previously identified from K imaging 
to be at z - 0.8 (Haines et al. 2001a; Haines 2001) as well as a 
z = 1 -334 quasar from the LQG. 

The ultimate aim of the near-infrared imaging is to iden- 
tify galaxies with each cluster candidate irrespective of star- 
formation history to M* + 2 through photometric redshift es- 
timates based on the VIJK S data. Here we only consider the 
early-type galaxies which allow the redshift and extent of the 
clusters to be determined through the red sequence method. 

5.1. Colour-Magnitude Relations 

Figure 8 shows the / - K s colour of the cluster red sequence 
as a function of redshift. The solid curves show the predicted 
colour evolution of early-type galaxies as modelled by stellar 
populations formed in an instantaneous burst at Zf = 3.5 (thick) 
and Zf = 2.0 (thin). The points correspond to the colour of the 
red sequence for several known clusters at z ~ 1 (Benftez et 
al. 1999; Rosati et al. 1999; Stanford et al. 2002; Tanaka et al. 
2000) and those from the surveys of SED98, A93 and N01 . It is 
apparent that / - K s photometry allows the redshift of clusters 
out to z - 1.3 to be determined efficiently, with red sequences 
expected at / - K s 3.4 and / - K s 4.1-4.5 for clusters at 
Z = 0.8 and z = 1 .2 respectively. 
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Fig. 8. The / - K s colours of the cluster red sequence as a function 
of redshift. The solid curves predict the colour evolution of early-type 
galaxies as modelled by stellar populations formed in an instantaneous 
burst at Zf = 3.5 (thick) and z/ = 2.0 (thin). The points correspond to 
the colour of the red sequence for several known clusters at z - 1 as 
well as clusters from the surveys of SED98, A93 and N01. 



We identify passively-evolving galaxies at z — 0.8 as hav- 
ing V - 1 > 2.3 and 3.0 < / - K s < 3.75. These criteria were 
chosen partly a posteriori to cover the observed VIK S colours 
of red galaxies in fields 2, 3 and 4, but a comparison with 
Figs. 2 and 8 shows that these criteria effectively select passive- 
evolving galaxies at 0.8 ± 0.2, the uncertainty corresponding to 
that expected for a redshift estimate for a galaxy at z ~ 0.8 from 
its VIJK S colour. 

Passively-evolving galaxies at z - 1 .2 are identified in the 
NTT fields as those having V - / > 1 .8 and 3.6 < / - K s < 4.8. 
The validity of these criteria can be confirmed, in a similar 
way as above, through comparison with the predicted colours 
of passively-evolving galaxies at 1 .2 ± 0.2 shown in Figs. 2 
and 8, but also with the results of the Gemini Deep Deep 
Survey (GDDS; Abraham et al. 2004). GDDS aimed at ob- 
taining a mass-selected sample of galaxies at 1 < z < 2, by 
preferentially selecting those near the / - K s versus / colour- 
magnitude track mapped out by passively evolving galaxies 
in this redshift interval, and succeeded to secure high-quality 
spectra of quiescent galaxies at 1 < z < 2. They identified 
twelve passively-evolving galaxies at 1 .2 < z < 1 .4, all of them 
having V - 1 > 1.8 and 3.5 </ - K s < 4.8. In contrast all of 
their z < 0.8 galaxies, excepting one, have I - K s < 3.6. The 
characterisation of a bluer V - I selection criterion for z - 1 .2 
galaxies with respect to those at z - 0.8 can be qualitatively un- 
derstood from Fig. 2. In effect, the change on the galaxy V - I 
colours due to the presence of star-formation, even in a small 
amount, is much larger for the z - 1.2 galaxies; this is appar- 
ent from the divergence at z - 1 .2 of the curves corresponding 
to the models describing passive and disk-dominated galax- 
ies at z - 1.2. Such behaviour is confirmed by the GDDS, as 
five of the twelve galaxies at 1 .2 < z < 1 .4 spectroscopically- 
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Fig. 9. Spatial distribution of galaxies in the four NTT fields. Each open circle indicates the location of a galaxy with the colours expected 
of early-type galaxies at the redshift of the cluster identified in each field, whose size indicates the galaxy's K s magnitude, the area being 
proportional to the observed flux. For field 1 this corresponds to z — 1.2, and the galaxies are selected to have V - I > 1.8 and 3.6 < / - K s < 4.8, 
whereas for fields 2, 3 and 4 this corresponds to z — 0.8 and the galaxies are selected to have V - I > 2.3 and 3.00 < / - K s < 3.75. In field 4 
passively-evolving galaxies at z - 1.2 are indicated by solid symbols. The nominal cluster regions used in Fig. 10 are shown by large circles. 
The location of each quasar is indicated by a star symbol. 



classified as passively-evolving have 1.8 < V - I < 2.3. We use 
a slightly redder minimum in our / - K s selection than the one 
resulting from the GDDS, in order to minimise the foreground 
contamination expected in our case from galaxies at z - 0.8. 

Figure 9 shows the spatial distribution of those galaxies in 
the four NTT fields identified as early-type galaxies at the red- 
shift of the cluster located in each field (z - 1.2 for field 1; 
z - 0.8 for fields 2, 3 and 4) from their photometry. Each galaxy 
is indicated by an open circle, whose size indicates the galaxy's 
K s magnitude, the area being proportional to the observed flux. 
The large circles indicate the nominal cluster regions used in 
Figs. 7 and 10. The location of each quasar is indicated by a 
star symbol. 

Figure 10 shows the / - K s against K s colour-magnitude 
diagrams for galaxies in the four NTT fields. To highlight the 
signal from the cluster red sequences in each of the fields, cir- 
cular regions around each of the nominal cluster centres are 
defined as shown in Fig. 9. For fields 1 and 4, one circular re- 
gion of radius 1' is defined, and for fields 2 and 3 where there 
is an apparent bimodal structure, 2 circular regions of radius 
(2~ 1/2 )' are defined. Galaxies within these circular regions are 
indicated by squares, with those having (not having) the colours 
of early-type galaxies at the cluster redshift indicated by solid 
(open) symbols. The remaining galaxies within the NTT fields 
are indicated by small crosses. 



Field 1 contains the z - 1 .2 cluster previously identified by 
Haines et al. (2001b), and this is confirmed by the excess of 
galaxies with I - K s ^ 4 in the circular region around the nom- 
inal cluster centre. This is as expected given that the circular 
region corresponds approximately to the 2.25 x 2.25 arcmin 2 
UFTI field used to obtain the K imaging described in Haines 
et al. (2001b). Although a tight red sequence is not apparent, 
28 K s < 20.5 galaxies within 1 arcmin of the cluster centre are 
classed as early-type galaxies at z - 1 .2 from their photometry, 
whereas only 2-4 would be expected. Figure 9 (top-left) shows 
that the clustering appears significantly extended beyond that 
identified by Haines et al. (2001b), with ~ 80 galaxies having 
colours consistent with passively-evolving galaxies at z - 1 .2 
across the NTT field, forming an amorphous, clumpy structure 
4-5 /^'Mpc in extent. 

In each of the three remaining fields, clear red sequences of 
20 galaxies are apparent at / - K s 3.4-3.5 in the colour- 
magnitude distributions of galaxies in the circular regions 
around the nominal cluster centres. Indeed the red sequences 
for the clusters in fields 3 and 4 appear almost identical, with 
that of field 2 marginally redder (Al - K s ~ 0.1), and com- 
parison with Fig. 8 indicates that all three clusters are at 
z = 0.8 + 0.1. The red sequence galaxies in fields 2 and 3 ap- 
pear bimodally distributed, with two compact groups of galax- 
ies located 3 and 5 arcmin apart respectively (Fig. 9). In con- 
trast the distribution of red sequence galaxies in field 4 appears 
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Fig. 10. The / - K s against K s colour-magnitude diagrams for galaxies in the four NTT fields. Galaxies in the circular regions corresponding 
to the clusters identified in each field are indicated by squares, with those having (not having) the colours expected of early-type galax- 
ies at the cluster redshift indicated by solid (open) symbols. For field 1 this corresponds to z — 1.2, and the galaxies are selected to have 
V - I > 1.8 and 3.6 < I - K s < 4.8, whereas for fields 2, 3 and 4 this corresponds to z - 0.8 and the galaxies are selected to have V - I > 2.3 
and 3.00 < I - K s < 3.75. The remaining galaxies in the fields are indicated by small crosses. 



unimodal, forming a single compact cluster with no obvious 
large-scale filamentary structure elsewhere in the field. 

In field 1, as well as the clustering at z ^ 1.2 manifested 
by the galaxies with / - K s - 4.1 (Fig. 10; top-left), there ap- 
pears a second distinct red sequence at / - K s ^ 3.4, identical 
in colour to those apparent in the remaining three fields, plac- 
ing them at the same redshift. In contrast to fields 2, 3 and 4, the 
galaxies with / - K s ^ 3.4 in field 1 appear evenly distributed 
across the whole NTT field. 

Although there is no further evidence of rich clusters at 
z - 1 .2 in the NTT data, in each of fields 2, 3 and 4, there are 
~ 10 K s < 19 galaxies with the VIJK S colours consistent with 
L*-class galaxies at z - 1 .2 with predominately old stellar pop- 
ulations. These do not appear clustered, except for a group of 
~ 6 galaxies with / - K s ^ 4.1 located near the z = 1.334 LQG 
quasar in field 4 (indicated by solid symbols in Fig. 9), suggest- 
ing that the latter is associated with this group. 

5.2. Star-Galaxy Separation 

The addition of near-infrared imaging allows stars and galax- 
ies to be distinguished on the basis of colour as well as mor- 
phology, with stars constrained to narrow loci in the J - K s or 
/ - K s against V - I colour-colour diagrams, appearing much 



bluer in the near-infrared than galaxies of comparable optical 
colours. It is thus possible to identify stars on the basis of their 
near-infrared colours, and to examine the efficiency of the mor- 
phological classification and estimating the level of stellar con- 
tamination in the red galaxy subset. We find that 10 of the 122 
sources in the three deepest NTT fields which are classified as 
stars on the basis of morphology have the near-infrared colours 
of galaxies; these are mostly faint (/ > 23) blue galaxies. Of the 
227 red galaxies in NTT fields 1-3 (field 4 is not used due to 
its insufficient depth), 17 are found to have the near-infrared 
colours of stars. By assuming that this density of stellar con- 
tamination is maintained over the whole BTC field, we would 
expect a total stellar contamination level of 250 ± 65 or 10%, 
where the quoted uncertainty assumes Poisson noise for the 17 
red galaxies identified as stars in the three NTT fields. As the 
image quality of the BTC and DLS images are similar in terms 
of depth and seeing levels, and lie at similar galactic latitudes 
(+53° for BTC; -50° for DLS), we would expect a similar level 
of stellar contamination in the DLS data. Therefore, although 
the number appears significant in comparison to the total ex- 
cess, the uncertainty in the relative levels of stellar contamina- 
tion between the two sets of data is just cr a 65 x V2 « 90, and 
so it is unlikely that more than 20% of the 563 galaxy excess 
(see section 4.2) is due to stellar contamination, particularly as 
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the magnitude and colours distributions of the red galaxy ex- 
cess and stellar populations are so different. 



6. DISCUSSION 

The realization that coherent massive structures > 100 /r'Mpc 
in size exist in any sufficiently large volume in the universe up 
to at least to z ~ 1 would provide a new major constraint for 
current and future theoretical models of the universe and its 
evolution. It is well known and widely accepted that structures 
50-100 /r'Mpc in diameter exist in the local universe, how- 
ever at higher redshifts it becomes much more difficult to detect 
similar structures and eventually we may be biased to detect 
mostly the largest ones. Very recently, two 100/i~'Mpc-scale 
structures at z ~ 0.3 were detected in the three-dimensional 
distribution of radio galaxies in a 25 deg 2 area (Brand et al. 
2003). In another 25 deg 2 region of the sky, previous investiga- 
tions have identified two 100 /r'Mpc-scale structures at z - 0.8 
and z - 1.2 (Clowes & Campusano 1991, 1994; Williger et al. 
2002), using counts of quasars and Mgn absorbers. It is not 
completely clear whether the volumes surveyed in these sur- 
veys are typical, or rather special, containing one or more rare 
fluctuations in the density distribution of the objects studied. 
In the case studied first by Clowes & Campusano(1991), al- 
though there are reasons to believe that an excess in the num- 
ber of quasars is marking a mass overdensity, the fact that it is 
also marked by an excess of Mg n absorbers is more directly 
indicative of an accompanying galaxy overdensity. In order 
to further characterize the distribution of the galaxy overden- 
sity implied by lOO/r'Mpc-scale structures in the Clowes & 
Campusano(1991) area, we have undertaken an ultra-deep op- 
tical imaging survey in a subarea of 40 x 35 arcmin 2 - denom- 
inated the BTC field -, to search for evidence of galaxy over- 
densities at z - 0.8 and z - 1 .2, either in the form of clusters, 
structures or a uniform distribution of galaxies at each redshift. 
Our analysis of the BTC data provided suggestive evidence for 
large-scale structure in this area of the sky, and allowed the 
identification of four cluster candidates all at z ~ 1. Next, we 
obtained targeted NIR imaging around these cluster candidates 
to ascertain the actual nature of the clusters and the large-scale 
structures, and to evaluate if they can be related with the quasar 
and Mg n absorber superstructures. Using both the optical and 
NIR data, we obtain results showing that, at least in the probed 
volumes, the two z ~ 1 100 h~ l Mpc-scale structures do indeed 
present an enhanced density of galaxies. As the BTC field is a 
representative subarea of the region where the superstructures 
lie, we consider our results to be evidence for the presence of 
sheet-like galaxy structures at z ~ 1 and an indication that they 
may be "Great Walf'-like structures at this early epoch of the 
universe. Although we describe the two z ~ 1 structures as 
"sheets", given the sparse spatial sampling of the superstruc- 
tures offered by the quasars and Mg n absorbers, and the lack 
of redshift information for the galaxies, we cannot accurately 
measure the topology of the superstructures. However, their ex- 
tent over the plane of the sky indicates that they are at least 
two-dimensional (as opposed to filamentary) structures. 



6. 1 . Selection of early-type galaxies at0.5<z<1.3 

To search for a possible large-scale structure of galaxies at 
z ~ 1 we have applied a variant of the cluster red sequence 
algorithm of Gladders & Yee (2000) to ultra-deep V and / 
imaging data of the BTC field. Instead of considering narrow 
redshift-slices over the redshift range of interest, a single high- 
redshift slice is used, containing all / < 23.5 galaxies redder in 
V - I than the cluster red sequence at z = 0.5. Simple galaxy 
evolution models predict that this selection of red galaxies ef- 
fectively identifies a population of passively-evolving early- 
type galaxies at 0.5 <z< 1.3, whose strong 4000A breaks are 
responsible for their red V - I colour. The validity of this ap- 
proach is demonstrated by application of the colour-magnitude 
criteria to galaxies in the EIS HDF-South field for which deep 
UBVRI (and for 40% of the field JHK S ) photometry exists. Of 
the 34 red galaxies identified: 32 were constrained at the l<x 
level to be at z > 0.6 from their photometric redshifts, the re- 
maining two constrained to z > 0.45; and all were best-fit by 
early-type galaxy models. 

6.2. A new large-scale structure of galaxies atz- 0.8 

Any large-scale galaxy structure at z ~ 1 underlying either of 
the quasar/Mg n superstructures at z - 0.8 or z - 1 .2 should 
manifest itself as an excess of red galaxies across a sub-area 
under study such as the whole BTC field. In effect, through 
comparison with a suitable control field (Flp22 from the Deep 
Lens Survey of Wittman et al. 2002) an excess of 563 red 
galaxies across the BTC field is identified. The magnitude and 
colour distributions of the excess red galaxies show unambigu- 
ously that the excess is largely due to a population of early-type 
galaxies at z - 0.8, and cannot be explained as either a calibra- 
tion error or due to stellar contamination. The colour distribu- 
tion of this excess shows a coherent peak at 2.7 < V - 1 < 3.1, 
the signature predicted for a population of early-type galax- 
ies at 0.5 <z< 1.3, and the magnitude distribution is well de- 
scribed by a Schechter function with a shallow-faint end slope 
(a = -0.3 ± 0.25) indicative of early-type galaxies (Madgwick 
et al. 2002). A comparison of m* with those of cluster popula- 
tions out to z - 1 (e.g. Nelson et al. 2001) places the excess at 
z = 0.83 ± 0.08. 

The presence and redshift of this z - 0.8 structure is con- 
firmed by the targeted near-infrared imaging of the four best 
galaxy cluster candidates that had been identified in the BTC 
field, resulting in the identification of three clusters at z = 
0.8 ±0.1, each having clear red sequences of ~ 20 galaxies 
with I - K s =! 3.4-3.5. Independent of the overall red galaxy 
excess, these clusters form a superstructure at z - 0.8 that is 
~ 20/T 1 Mpc across, which extends fully across the BTC field. 

There is no obvious apparent filamentary structure connect- 
ing the z =2 0.8 clusters in the spatial distribution of red galax- 
ies (Fig. 9). However as the clusters make up only 20% of the 
overall red galaxy excess, the most important contributor to the 
remainder is likely to be due to a connecting filamentary or 
sheet-like large-scale structure. Such a filamentary structure is 
difficult to detect from the red galaxy distribution, as it is likely 
to be diffuse, and the fraction of early-type galaxies will be 
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lower than in the clusters, reducing the signal still further. Only 
further near-infrared imaging or multi-object spectroscopy over 
a wide area is likely to detect such filamentary structures. There 
is however evidence of such a diffuse, filamentary or sheet-like 
structure in field 1, where despite there being no clusters at 
z - 0.8, a red sequence at / - K ^ 3.4 is apparent, made up of 
galaxies evenly distributed across the NTT field. This suggests 
that the overall excess of red galaxies is due to two main com- 
ponents, on the one hand, the three z — 0.8 clusters, and on the 
other hand, a diffuse sheet-like structure extending across the 
whole BTC field. 

6.3. A new large-scale structure of galaxies at z- 1.2 

Previous TT-band imaging of a 2.25 x 2.25 arcmin 2 field centred 
on the z = 1 -233 LQG quasar indicated the presence of associ- 
ated clustering in the form of red sequences of 15-18 galax- 
ies with the colours expected for a population of passively- 
evolving massive ellipticals at the quasar redshift (7 - K ^ 4.3, 
V - K =s 6.9). This field makes up part of NTT field 1 and as 
such has been reobserved with higher quality NIR data, with 
a field of view five times larger than before, having a mag- 
nitude greater depth in K, along with new /-band imaging. 
The previously observed red sequence is found again, with 28 
K s < 20.5 galaxies within 1 arcmin of the cluster centre, al- 
though now at I - K s ^4.1. This appears due to a systematic 
offset between the two data sets, and given the extensive and 
well-behaved photometric calibrations of the NTT data, we 
consider these to be correct. It had been suggested in Haines 
et al. (2001b) that the clustering extended beyond the K im- 
age, by consideration of the spatial distribution of optically 
red galaxies (Fig. 7 of Haines et al. 2001b), with a filament 
predicted to extend 3 arcmin to the east of the quasar, and a 
compact group 3 arcmin to the south. These are covered by the 
larger NTT field, and are apparent in Fig. 9(top-left), along 
with another group 2.5 arcmin to the west of the quasar. The 
overall structure extends across the NTT field, corresponding 
to 3-4/z~'Mpc at z — 1.2, and appears clumpy and filamen- 
tary. This appears to be a very massive system, considering its 
redshift, with ~ 80 galaxies (40 with K s < 19) having colours 
consistent with being passively-evolving massive ellipticals at 
Z — 1.2. There are also comparable numbers of galaxies having 
colours consistent with being star-forming galaxies at the same 
redshift (I - K s > 3.6, V - I < 1.8) as seen also for other z > 1 
clusters (Tanaka et al. 2000) and NIR surveys (McCarthy et al. 
2001). 

Although there is no further evidence of rich clusters at 
Z - 1 .2 in the BTC field, there is a group of ~ 6 galaxies with 
/ - K s =s 4.1 located near the z = 1.334 LQG quasar in field 4, 
suggesting that the latter is associated with this group. Also 
~ 10 galaxies with I - K s ^4.1 are apparent in each of the 
NTT fields 2, 3 and 4, which could indicate a diffuse structure 
at z - 1 -2 across the BTC field. 

It is difficult to ascertain the full extent of the z - 1 .2 large- 
scale galaxy structure in the regions where only optical data ex- 
ists, and to compare the structure with that of the z - 0.8 struc- 
ture. There are a number of effects which combine to signifi- 



cantly reduce the efficiency of the red galaxy selection criteria 
in detecting structures at z = 1.2 as opposed to z = 0.8. Firstly 
the galaxies are 1-1.5 mag fainter due to the greater distance 
and the effects of k-corrections. Secondly the V - 1 colour 
samples shorter rest-frame wavelengths (from 3000^-500A at 
Z = 0.8 to 2500-3600A at z = 1.2), so the same amount of re- 
cent star-formation has approximately double the effect on the 
V - 1 colour of a galaxy at z = 1 .2 as it would in the same 
galaxy at z = 0.8. The galaxies themselves are younger (by 
1.2/T 1 Gyr) and so are intrinsically bluer. Finally as the clus- 
ters themselves are younger and less dynamically evolved, their 
inhibitory effect on star-formation in their member galaxies is 
reduced, as is evident from the observations of the Butcher- 
Oemler (1984) effect. To date the redshift limit to which cluster 
red sequences have been found is z ~ 1.3. 

Many of these effects can be partially counteracted by the 
addition of NIR photometry, which is less affected by star- 
formation, allowing a greater fraction of the early-type galax- 
ies to be detected. Of the ~ 30 galaxies in fields 2, 3 and 4 
with I - K s ~ 4.1 (and hence candidate early-type galaxies at 
Z - 1 .2), only ~ 50% would be classed as red galaxies, either 
appearing too blue in V - I (either due to some star-formation 
or photometric uncertainties as the limiting magnitude in V is 
reached), or have / > 23.5 and are thus too faint. Hence ap- 
proximately half the galaxies likely to be early-type galaxies at 
Z - 1 .2 are missed by the optical selection criteria, and it is only 
the very rare high-density peaks corresponding to rich clus- 
ters that become apparent at z - 1 .2 in the red galaxy spatial 
distribution. In contrast it is only in the regions with the NIR 
data that the sheet-like large-scale structure associated with the 
Clowes-Campusano LQG becomes apparent. This loss of effi- 
ciency in identifying early-type galaxies at z > 1 from optical 
data alone is apparent from the results of the Las Campanas 
Infrared Survey (McCarthy et al. 2001) where galaxies identi- 
fied as being the progenitors of early-type galaxies at 1 < z < 2 
from having I - H > 3 span a range in V - I colour of more 
than 3 mag, with around half of them having V - I < 2. 



6.4. Large-scale structures atz < 0.3 

To date the distribution of galaxies on 100/r'Mpc-scales has 
been followed out to z ~ 0.3 (e.g. Geller & Huchra 1989; 
Doroshkevich et al. 2000), which is also the upper limit of 
the major galaxy redshift surveys such as the 2dF GRS survey 
(Colless et al. 2001) and the Sloan Digital Sky Survey (SDSS; 
York et al. 2000). In consequence studies into the relation be- 
tween quasars and the large-scale galaxy distribution have been 
limited to z < 0.3. Sochting et al. (2002, 2003), although based 
on 2-colour photographic sky survey data, have presented not 
only evidence of large-scale structure at z - 0.3 delineated by 
galaxy clusters but also that quasars occur preferentially close 
to these structures. Two 100/z _1 Mpc structures of radio galax- 
ies have been recently discovered at z ~ 0.3 in a 25 deg 2 area 
of the sky (Brand et al. 2003), which are estimated to be highly 
improbable under the quasi-linear structure formation theory 
and thus need to invoke an increased bias on the large scales. 
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6.5. Large Quasar Groups as tracers of large-scale 
structure at high-redshifts 

The detection here of large-scale structures of galaxies that de- 
lineate part of known quasar/Mg n superstructures at z - 1, is 
similar to that found by Tanaka et al. (2000, 2001) for the 
Crampton et al. (1989) group of 23 quasars at z - 1.1. In a 
study comparable to that of Haines et al. (2001b), Tanaka 
et al. (2000) obtained deep R-, I- and TT-band imaging of 
a 5 x 3.2arcmin 2 region centred on the z = 1.086 radio-loud 
quasar 1335.8+2834 from the Crampton et al. LQG. A rich 
cluster of galaxies with the colours predicted of passively- 
evolving galaxies at z ^ 1 . 1 (/ - K ^ 4.0, R - K ^ 5 .3) is iden- 
tified, and, as is found by Haines et al. (2001b), the quasar is 
located on the cluster periphery in a star-forming region, in- 
dicated by a band of emission-line galaxies (Hutchings et al. 
1995), and also there is evidence of further filamentary struc- 
tures forming a large-scale structure 4-5 h~ 1 Mpc in extent. 

In a wide-area (48 x 9 arcmin 2 ) imaging survey in both R- 
and /-bands toward a region containing five quasars from the 
Crampton et al. (1989) LQG, Tanaka et al. (2001) detect sig- 
nificant clustering of faint, red galaxies with 21 < / < 23.5 and 
1 .2 < R - I < 1.6, i.e. galaxies with the colours and magni- 
tudes expected of early-type galaxies at z - 1.1. These galax- 
ies are concentrated in 4-5 clusters forming a linear struc- 
ture of extent ~ 10/z^'Mpc that is traced by members of 
this large group of quasars, although only the one radio- 
loud quasar of Tanaka et al. (2000) appears directly as- 
sociated with any of the rich clusters. Conversely, in an 
optical spectroscopic follow-up of six X-ray sources de- 
tected in a 20 x 20 arcmin 2 field, known to contain two rich 
galaxy clusters at z ~ 1.3 (RXJ0848.6+4453 at z = 1.273 and 
RXJ0848.9+4452 at z = 1.261; Stanford et al. 1997; Rosati et 
al. 1999), three of the sources were identified as quasars at 
z = 1.260, z = 1.286 and z = 1.260 (Ohta et al. 2003). 

There are two main factors which should affect how quasars 
trace mass: the requirement of some disturbance to the host 
galaxy to push gas onto the central nucleus and activate the 
quasar, either due to a merging event, which is likely to be the 
dominant mechanism at z - 1, or during the formation of the 
galaxy (Haehnelt & Rees 1993); and the requirement of a mas- 
sive black-hole to fuel the quasar, which given the strong cor- 
relation between the masses of the black hole and the bulge of 
the host galaxy (Merritt & Ferrarese 2001), implies a massive 
host galaxy. 

The effect of quasars being located in merging/forming 
galaxies on how they trace mass is likely to be complex and 
redshift dependent. However on large-scales the effect should 
be small, and the most important consequence should be the 
avoidance by quasars of the high-density cluster centres at 
least to z ~ 1 • This is understandable in the framework of both 
galaxy merger and galaxy formation quasar triggering mecha- 
nisms: the encounter velocities of galaxies in the centre of clus- 
ters are much greater than the internal velocity distributions, 
and so galaxy mergers become much less efficient at triggering 
nuclear activity (Aarseth & Fall 1980); and the cluster cores are 
filled with shock-heated virialised gas that does not easily cool 
and collapse (Blanton et al. 1999), inhibiting both the forma- 



tion of stars and galaxies (Blanton et al. 2000), and hence in- 
hibiting quasar formation. This effect has been observed in the 
form of the preferential location of quasars on the peripheries 
of clusters (Sanchez & Gonzalez -Serrano 1999, 2002; Haines 
et al. 2001b; Sochting et al. 2002; Tanaka et al. 2000). 

The observation that luminous (My < -23.5) quasars (at 
least to z ~ 1) are located in massive ellipticals (e.g. McClure et 
al. 1999; Kukula et al. 2001), and the requirement of a massive 
host galaxy, indicates that on large-scales at least ( > 5 /z^Mpc) 
quasars should trace mass in the same way as their host galax- 
ies, a prediction that appears confirmed by comparison of 
the quasar and galaxy power-spectra (Hoyle et al. 2002). We 
should thus expect quasar superstructures to trace the same 
mass overdensities as the large-scale structures mapped by 
galaxies, as is suggested by comparison of the quasar and 
galaxy spatial distributions in the local universe (Longo 1991). 

The results presented here and elsewhere (Tanaka et al. 
2001; Ohta et al. 2003) indicate that super-large scale struc- 
tures (~ 100/r'Mpc) in the form of galaxies, comparable to 
the "Great Wall" in the local universe, are apparent at z ~ 1, 
and that they can be discovered and traced well by quasar 
superstructures. This suggests that large-scale quasar surveys, 
such as the 2dF QSO survey, could be used to map the evo- 
lution of the super-large scale structure to redshifts well be- 
yond the scope of current galaxy redshift surveys, reaching 
z - 1 or beyond. Statistically significant structures on scales 
> lOO/z^Mpc have been identified in the 2dF QSO survey 
(Miller et al. 2004), appearing rarely, but throughout the sur- 
vey, and over the whole redshift range covered (0.5 < z < 2.2). 

6.6. Evolution of Large-scale Structures 

The finding of large-scale structures at z ~ 1 can be under- 
stood in the context of the Zel'dovich non-linear theory, as the 
perturbations in the initial density field collapse first to form 
pancakes - the sheet-like structures of the LSS - before then 
collapsing on the remaining two axes to form filaments and 
clusters. As clusters and filaments are apparent at z ~ 1, the 
sheet-like structures which form before these, must also be sub- 
stantially in place. This is confirmed by simulations, in which 
Doroshkevich et al. (1999) find that for the present day wall- 
like structures where the overdensities of matter are a factor 
5-10 over the mean, half the mass was already in structures 
with the same level of overdensities at z ~ 1 . This also can be 
considered through the simple argument that the characteris- 
tic distance scales of large-scale structures (50-100 /z _1 Mpc) 
are much greater than the distances galaxies are likely to move 
over the time-scales involved, (a galaxy with a peculiar velocity 
of 100 hT x kms -1 moves 1 h~ l Mpc over lOGyr). Hence, those 
galaxies that make up the large-scale structures apparent today 
and which were formed at z > 1, will still constitute the same 
large-scale structures at z ~ 1 . 

The observation of two sheet-like large-scale structures 
at z ~ 1 in the same field, in conjunction with the results of 
Tanaka et al. (2001) and Ohta et al. (2003), implies that such 
structures although rare, are not exceptional. The observed 
early formation and relatively slow evolution of the abundances 
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of massive structures, clusters in particular, is one of the funda- 
mental reasons for favouring the current ACDM models (e.g. 
Bahcall & Bode 2002). The Hubble Volume numerical simula- 
tions of Evrard et al. (2002) find a large cluster at z = 1 .04 in 
a ACDM model which has a mass twice that of Coma, making 
it the largest cluster in their positive octant (PO) survey which 
covered an eighth of the sky to z = 1 .46. In numerical simula- 
tions Doroshkevich et al. (1999) find that sheet-like large-scale 
structures with typical maximal extensions of ~ 30-50 /z~'Mpc 
incorporate ~ 40% of matter at z = 0, of which ~ 60% is in 
place by z ~ 1 . 

7. CONCLUSIONS 

The results presented here support the presence of two 
sheet-like large-scale galaxy structures extending across the 
20-30 /r'Mpc scales explored by our optical data: one "sheet" 
at z — 0.8 with three associated clusters; and a second "sheet" 
at z — 1 .2 with one embedded cluster. Therefore, our results 
confirm that, at least in the volumes probed, the two z ~ 1 
structures mark out volumes with an enhanced density of 
galaxies. As the studied subarea is quite representative of 
the whole region where the superstructures lie, we consider 
these results to be evidence for the presence of large sheets 
of galaxies at z ~ 1, and a clear indication that such sheets 
may extend over 100/r'Mpc-scales. Secondly, these results 
show that large quasar groups (LQGs) are reliable indica- 
tors of galaxy excesses at ~100Mpc-scales. In particular, the 
Clowes & Campusano (1991) LQG, the largest structure at 
high redshifts yet found, is then indicative of a huge struc- 
ture of galaxies (and thus of coherent mass) over a volume of 
~ 100 x 200 x 200 /T 3 Mpc 3 . Although the observational deter- 
mination of the spatial frequency of such fluctuations in the 
distribution of galaxies at z ~ 1 will require future generations 
of galaxy surveys, the Clowes-Campusano LQG provides an 
extraordinary laboratory for the study of the relation of struc- 
tures at different scales; the interplay between radio-galaxies, 
quasars and gas; and the evaluation of the bias factor for dif- 
ferent constituents and environments inside the superstructure. 
Thirdly, the handful of LQGs that have been discovered so far 
at 0.2 < z < 2.0 (see Graham et al. 1995), suggest that super- 
structures of this size are indeed rare, corresponding to unusu- 
ally high peaks in the primordial density field. Brand et al. 's 
(2003) discovery of two lOO/r'Mpc-scale structures of radio 
galaxies at z ~ 0.3 may also correspond to related evolution of 
similarly rare peaks. 

The detections of two lOO/z'Mpc-scale structures in sur- 
veyed volumes such as that of Clowes & Campusano and Brand 
et al. studies, was shown by the latter to be incompatible with 
DM models predictions if these volumes are typical. Brand et 
al. explained the observations as a consequence of an increased 
"bias factor" and concluded that structures 100/r'Mpc in size 
at z ~ 0.3 do not present a challenge to the standard inflationary 
ACDM model. Using a semi-analytic model which includes 
both galaxy and quasar formation, in which the quasar activ- 
ity is fuelled by major galaxy mergers, Enoki et al. (2003) 
find that quasars at z ~ 1 .2 are biased with respect to the un- 
derlying mass distribution, with b ~ 2, and are found in a va- 



riety of environments from small galaxy groups to rich clus- 
ters, as is observed (e.g. Wold et al. 2000, 2001). Equally, 
the massive galaxies which we are observing here, and which 
host the quasars in the LQGs, should have been twice as "bi- 
ased" at z = 1 as they are at z = 0, having b ~ 2 (Blanton et 
al. 2000). These bias factors allow the finding of such super- 
structures of quasars and galaxies to be reconciled with ACDM 
models. Further theoretical work on this topic, together with 
specific predictions for the size distribution of the super-large 
structures as a function of redshift, is highly desirable. By 
assuming a quasar "bias" redshift dependence similar to that 
described above Miller et al. (2004) find that the significant 
> 100/z'Mpc-scale fluctuations in the quasar distribution ob- 
served in the 2dF QSO survey are in good agreement with pre- 
dictions of a standard ACDM model. 

As the ability of infra-red cameras to cover large areas of 
sky (tens of deg 2 ) is expected to improve greatly over the com- 
ing years (e.g. UKIRT-WFCAM, VISTA), so will the possi- 
bility of doing multi-colour studies of the evolution and the 
bias function of super-large scale structure to z ~ 1 and beyond. 
Also the newest generation of wide field spectrographs hav- 
ing large multiplexing (100-4000), either built (ESO-VIMOS 
and NIRMOS, and Magellan-IMACS and Keck-DEIMOS) or 
planned (Gemini-Optical MOS), allow the planning of redshift 
surveys for the study of specific populations of galaxies belong- 
ing to the largest structures yet known. 
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